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A burnup calculation was performed to analyze the Aprés ORIENT process, which aims to create highly- 
valuable elements from fission products separated from spent nuclear fuels. The basic idea is to use nuclear 
transmutation induced by a neutron capture reaction followed by a B~ decay, thus changing the atomic number 
Z of a target element in fission products by 1 unit. LWR (PWR) and FBR (MONJU) were considered as the 
transmutation devices. High rates of creation were obtained in some cases of platinum group metals (4,Ru 
by FBR, 4sPd by LWR) and rare earth (¢4Gd by LWR, 66Dy by FBR). Therefore, systems based on LWR and 
FBR have their own advantages depending on target elements. Furthermore, it was found that creation rates of 
even Z’ (= Z + 1) elements from odd Z ones were higher than the opposite cases. This creation rate of an 
element was interpreted in terms of “average 1-group neutron capture cross section of the corresponding target 
element (0c zy" defined in this work. General trends of the creation rate of an even (odd) Z’ element from the 
corresponding odd (even) Z one were found to be proportional to the 0.78th (0.63th) power of (o.z), however 
with noticeable dispersion. The difference in the powers in the above analysis was explained by the difference 
in the number of stable isotopes caused by the even-odd effect of Z. 
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I. INTRODUCTION 


Recently, demand for Pd, Nd, Dy, and other rare metals is 
increased because they are used for many rapidly-progressing 
technologies such as clean energy devices. As a result, in 
Japan, which relies on the import of them from other coun- 
tries, supply risks of them have increased. Thus, a stable sup- 
ply of these rare elements is strongly desired. 

It is known that platinum group metals (PGMs), rare earth 
(RE), and other useful elements are contained in spent nucle- 
ar fuels (SNF) of nuclear reactors as fission products (FPs). 
They are defined nuclear rare metals (NRMs) in the concept 
of the Advanced ORIENT cycle [1], shown in the upper part 
of Fig. 1, in which these materials are actively retrieved from 
SNF by chemical partitioning to use as precious resources. 
Most of these elements, however, have very high or long- 
life radioactivity which prevents them to be considered as re- 
sources. 

The nuclear transmutation process has a possibility of de- 
creasing the amount of radioactive nuclides, as well as to 
change an element to another one, eventually leading to a sta- 
ble isotope by such a process: 


99Te(n, y) Te E °Ru(stable). (1) 


The Aprés ORIENT research program [2] shown in the low- 
er part of Fig. 1 was initiated in FY2011 aiming at creating 
stable, highly-valuable elements (i.e. secondary NRMs) by 
nuclear transmutation from FPs contained in SNF of nucle- 
ar reactors. The present investigation was carried out as a 
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Fig. 1. (Color online) Scheme of Advanced ORIENT cycle concept 
and the Après ORIENT research program. 


fundamental study to embody the Après ORIENT research 
program to clarify characteristics of creation of elements by 
nuclear transmutation. 


Il. CALCULATION METHOD 


A burnup calculation was performed to analyze the nuclear 
transmutation rate from FPs. In this study, nuclear reactors 
were assumed as sources of intense neutron flux and neutron 
capture reactions by reactor neutrons were considered as a 
method of creation of NRMs by nuclear transmutation. Af- 
ter a nucleus captures a neutron, it may eventually undergo 
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B- decay and, therefore, is transformed into another element 
having one proton more than before. Hence, stable highly- 
valuable rare metal elements (i.e. secondary NRMs) will be 
created by (n, y) reaction, succeeded by B~ decay of FPs 
which are generally characterized by high radioactivity. Such 
a process can be written as 


AFP(n, y)^tIFP £ + 4+INRM(stable). (2) 


The burnup calculation was carried out with the ORLIBJ40 
package [4], which is a combination of ORIGEN2.2 code [5] 
as burnup calculation code and 1-group cross section library 
based on JENDL-4.0 [6]. In ORLIBJ40, the irradiation ge- 
ometry is assumed to be homogeneous, which may not be 
realistic. This point will be improved in our future work by 
employing 3-dimensional Monte-Carlo simulators. 


ORLIBJ40 
Calculation Code : ORIGEN2.2 
Cross-section Library : JENDL-4.0 


LWR PWR : Running (1125 days) 
(35U4.7% , 45000MWd/tHM) 


Cooling (5 years) 


392 ~ goTM Reprocessing 


Separating Elements from Fission Products 
Re-loading One Element with Fuel 


Capture cross-section is large. Neutron Flux is large. FBR 


PWR : Running (1125 days) MONJU : Running (1125 days) 
) | (Neutron Flux : 2.70 X 1015/cm?/s) 


(Neutron Flux : 3.29 X 10!4/cm?2/s 
Library : PWR47J40 = Library : MONJMXRDJ40 
Cooling (5 years) 


Separating Created Element from Target 


Fig. 2. (Color online) Flowchart of computation for the creation of 
secondary NRMs in Aprés ORIENT scheme. 


The flowchart of computation adopted in this study is 
shown in Fig. 2. It is assumed to irradiate a target element 
by neutrons in LWR such as PWR, FBR such as MONJU (a 
sodium cooled prototype fast breeder reactor in Japan). The 
computation scheme consists of the following 7 steps: 


STEP 1: A burnup calculation is performed with the condi- 
tion that a PWR is operated for 1125 days to reach 
45 000 MWd/tHM. The fuel is assumed to be fresh 
with an enrichment of 7°°U by 4.7%. A 1-group 
cross section library PWR47J40.LIB [7] is used in 
this calculation. 


STEP 2: A decay calculation of the spent nuclear fuels during 
a cooling period of 5 years is carried out. Yields 
of elements contained in FPs after the cooling are 
shown in Fig. 3. 


STEP 3: Amount and isotopic composition of each elemen- 
t are analyzed as an initial condition of the target 
zFP element for neutron irradiation in the follow- 
ing steps. 
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TABLE 1. Isotope composition of the nuclear fuel 


Mass number LWR (PWR) FBR (MONJU) 
(%) (%) 
U 235 4.7 0.24 
238 95.3 80.76 
Pu 238 — 0.57 
239 — 10.07 
240 = 4.75 
241 — 2.28 
242 — 1.33 


TABLE 2. Isotope composition of reloaded Ru and Rh as examples 
of reloaded a 7FP (*metastable) 


a4Ru asRh 
Mass Weight Half-life Mass Weight Half-life 
number (g) number (g) 
98 1.20E-07 Stable 101 3.54E-10 3.3 years 
99 2.35E-02 Stable 102 3.22E-06 207 days 
100 1.42E+02 Stable *102m 5.51E-08 3.742 years 
101 1.05E+03 Stable 103 6.13E+02 Stable 
102 1.07E+03 Stable 106 6.17E-06 30.07 sec 
103 5.41E-13 39.26 days 
104 7.21E+02 Stable 


106 6.59E+00 371.8 days 


STEP 4: Input data of each target ¿FP element is prepared 
for a burnup calculation of a 7FP assumed to have 
been obtained by perfect mutual element separation 
(amount and isotopic composition of which is deter- 
mined in STEP 3). Then, the 7FP is also assumed 
to have been reloaded in a core region of a PWR 
or a blanket region of MONJU at the time of the 
mutual element separation. Isotopic composition of 
the nuclear fuel of each reactor core is shown in Ta- 
ble 1. Isotopic composition of reloaded Ru and Rh 
are shown in Table 2 as examples of reloaded a FP. 


STEP 5: A burnup calculation is performed corresponding to 
operation of each reactor for 1125 days with the fol- 
lowing conditions: 


PWR; neutron flux: 3.29 x 10!4/(cm?s), with a 1- 
group cross section library as PWR47J40.LIB. 


MONJU; neutron flux: 2.70 x 10'°/(cm?s), with 
a l-group cross section library as MONJMXRD- 
J40.LIB [7]. 


STEP 6: A decay calculation is carried out corresponding to 
a cooling period of 5 years. 


STEP 7: Analysis of the amount and isotopic composition of 
the 7,,NRM element created from the FP is per- 
formed. 
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Zn Ga Ge As Se Br Kr Rb Sr Y Zr NbMoTc Ru Rh Pd Ag Cd In Sn Sb Te | Xe Cs Ba La Ce Pr NdPmSm Eu Gd Tb Dy Ho Er Tm 


(Color online) A bar chart showing calculated results on the amounts of FP elements in a PWR spent fuel from 1 t of fresh nuclear fuel 


(U: 4.7%) after an operation time of 1125 days (which leads to 45 000 MWd/tHM) followed by a cooling of 5 years. 
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Fig. 4. (Color online) Comparison of creation rates of each 7,,NRM element. 


RESULTS AND DISCUSSION 


A. Comparison of the “creation rate” of a 7,,;,NRM element 


In order to discuss transmutation efficiency quantitatively, 
“creation rate” of a 7,, NRM element was defined as 


( 


Creation rate 
ofz +ıNRM 


Amount of 


( created z ,,NRM ) H 


) (%/year] = 


( 


( 


Amount of 


re-roaded Te) le] 


100 


Irradiation 
time 


) [year] 


(3) 


Actually, this quantity indicates a yearly-average of the trans- 
mutation rate for an irradiation period of 1125 days (3.08 
years). Then, comparison of creation rates of each 7,,NRM 
element calculated by Eq. (3), as shown in Fig. 4. The cre- 
ation rate of ,;Rh represents the calculated result for ,;Rh 
created from transmutation of reloaded ,,Ru contained in F- 
Ps. It turned out that high creation rates appeared in valu- 
able elements such as PGMs (,,Ru: 8.03%/year (FBR), 4¢Pd: 
20.6%/year (LWR)) and middle or heavy RE elements (for 
example, ¢4Gd: 23.4%/year (LWR), 6¢Dy: 17.1%/year (F- 
BR)). Overall, there is a tendency in many elements having 
small creation rates (less than 1%/year) such as 3,Y, that the 
creation rate by FBR is higher than that by LWR, depend- 
ing on the difference in the amount of neutron flux. On the 
contrary, some elements with high creation rates (more than 
10%/year) such as 5Sn, show the opposite results regardless 
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Fig. 5. (Color online) Comparison of average 1-group neutron capture cross sections of each target zFP element. 


of the difference in the neutron flux. Therefore, it can be 
concluded that both transmutation systems have certain ad- 
vantages and disadvantages and it should be carefully chosen 
which system to use to create a certain 7,,NRM. Inciden- 
tally, if proper moderators were applied in the transmutation 
region of FBR, higher creation rates could be obtained for any 
zı NRM element. 

It is also clear in Fig. 4 that there is a difference in these 
creation rates according to even-oddness of atomic numbers 
Z' (= Z +1). In addition, it is shown that these elements 
with even Z’ have higher creation rates than those with odd 
Z'. For example, the creation rate of ,.Pd from ,;Rh is higher 
than those of both sides, that is, ,;Rh from ,,Ru and ,,Ag 
from ,,Pd. It is generally known that nuclei with an even 
number of protons are more stable than those with an odd 
number of protons. Therefore, the difference in each creation 
rate of 7,,NRM is considered to be based on the difference 
in the stability of the nucleus due to even-odd effects. 


B. Introduction of “average 1-group neutron capture cross 
section” 


The transmutation rate of a nuclide can be approximately 
represented to be directly proportional to the neutron capture 
cross section of the nuclide and neutron flux, such as 


transmutation rate capture cross section 
of a nuclide Z; of the nuclide ocz, (4) 


x (neutron flux ¢). 


Naively, the creation rate ofa 7,,NRM element might be rep- 
resented in the same way. However, the cross section of a 
zFP element does not exist because a cross section is a phys- 
ical quantity associated with a nucleus but not an element. In 
the present study, the “average 1-group neutron capture cross 
section” of a ¿FP element was introduced, and was defined 


as 


= > Ocz; MZ, 
(Jez) T Sma (5) 


(ocz): average l-group neutron capture cross section of 
element Z; 


Ocz,: 1-group neutron capture cross section of isotope 7 
of element Z; 


mz,: weight of isotope 2 of element Z in FP. 


Figure 5 shows average 1-group neutron capture cross sec- 
tions of each target 7FP element. It turns out that the struc- 
ture of these average 1-group neutron capture cross sections 
of each target FP element (Fig. 5) has a clear correlation 
with that of those creation rates of each 7,,NRM elemen- 
t (Fig. 4). In addition, it is noticed that there are two valleys 
minimizing cross sections (around Sr and Ce). This structure 
is considered to be formed by the presence of N = 50 and 82 
“neutron magic numbers” in this mass region, as indicated in 
Table 3. 


TABLE 3. Nuclides with neutron magic numbers in target 7FP ele- 
ments reloaded (*radioactive nuclide) 


Neutron magic number N =50 N = 82 
Kr 86 Xe 136 
Rb 87 Cs 137* 

The nuclide which exists in target Sr 88 Ba 138 
zFP elements reloaded Y 89 La 139 
Zr 90 Ce 140 

Nb 91* Pr 141 

Mo 92 Nd 142 
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Fig. 6. (Color online) Correlation between creation rates of each z,;NRM element and average 1-group neutron capture cross sections of 
these corresponding target zFP elements. (a) exhibits results of transmutation from odd Z target FPs to even Z’ NRMs. On the other hand, 


(b) exhibits from even Z target FPs to odd Z’ NRMs. 


C. Relation of creation rates of 7,,;,NRM elements with 
average 1-group neutron capture cross sections of target ZFP 
elements 


Figure 6(a) shows creation rates of even Z’ (= Z + 1) N- 
RMs as a function of average 1-group neutron capture cross 
sections of target FPs having an odd Z. Fig. 6(b) shows sim- 
ilar data for even Z target FPs leading to odd Z’ NMRs. Be- 
cause FBR (MONJU) had a larger neutron flux than LWR (P- 
WR), creation rates by FBR were higher than those by LWR. 
In addition, the larger the average 1-group neutron capture 
cross section of each target FP element was, the higher the 
creation rate of each 7,,NRM element was. Thus, a clear cor- 
relation was established between creation rates and average 
1-group neutron capture cross sections. In the case of the cre- 
ation of NRMs with even Z’, the creation rate of a 7,,NRM 
element was approximately proportional to the 0.78th power 
of the average 1-group neutron capture cross section of the 
target ZFP element, as 


creation rate of an element (o \0.78 (6) 
with even atomic number a i 


In addition, deviation from this trend is rather small. On the 
other hand, in the creation of NRMs with an odd Z’, the cre- 
ation rate of a 7,,NRM element was roughly proportional to 
the 0.63th power of the average 1-group neutron capture cross 
section of the corresponding target zFP element, as 


creation rate of an element (ez), (7) 
with odd atomic number 


However, this case is characterized with a large dispersion 
from the average trend. Reasons for the difference in the pow- 
ers and dispersion will be discussed in the following section. 


TABLE 4. The difference in the created 7,,NRM element by the 
even-odd effect of atomic number 


Atomic number of created 7,,NRM 


The number 
of stable isotopes 


Even [ex. Pd] 
lor2 
(It is easy to 


Odd [ex. Rh] 
Many 
(It is hard to 


of Target change elements.) change elements.) 
zFP [Rh:1] [Ru:7] 

The number Many lor2 

of stable isotopes (Created elements (Created elements 
of Created are stable.) are unstable.) 
zı NRM [Pd:6] [Rh:1] 

The dispersion 

from the average Small Large 


trend 


D. Classification of the difference in the even-odd effect of 
atomic number 


Table 4 summarizes the difference the even-odd effect has 
an atomic number Z’ (= Z + 1) from created NRMs. In the 
case of creation of NRMs with even atomic number Z’ (ex. 
aged), target FPs (,,Rh) have an odd atomic number Z, and 
they have only 1 or 2 stable isotopes (,;Rh: 1). That is, if 
an isotope (4;Rh) captures a neutron, it can decay easily, as 
in Fig. 7. Additionally, the 7,,NRM element (,,Pd), created 
by 6~ decay of an isotope of the target ,FP element cap- 
turing a neutron, has so many stable isotopes (,,Pd: 6) that 
it hardly causes additional transmutation to another elemen- 
t (ex. 4,Pd— ,,Ag). Fig. 8 shows the time variation of the 
weight of isotopes of target ,,Rh in the case of irradiation by 
LWR. It turns out that because a majority of isotopes of tar- 
get ,.Rh are almost all !°°Rh, !°'Rh is mainly transmuted and 
decreased, as 


103Rh(n, y)04Rh s 14Pa(stable). (8) 
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Fig. 7. (Color online) Main burnup chain from target 45Rh to created 46Pd. 
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Fig. 9. (Color online) Time variation of the weight of isotopes of 4¢Pd created from target 4;Rh in the case of irradiation by LWR. 
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Fig. 10. (Color online) Main burnup chain from target ,,Ru to created 45Rh. 
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Fig. 11. (Color online) Time variation of the weight of isotopes of target ,,Ru in the case of irradiation by LWR. 
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Fig. 12. (Color online) Time variation of the weight of isotopes of ,;Rh created from target 44Ru in the case of irradiation by LWR. 


In addition, Fig. 9 shows the time variation of the weight of 
isotopes of ,,Pd created from target ,,Rh in the case of irradi- 
ation by LWR. It turns out that many isotopes of ,,.Pd, which 
are heavier than 1%Pd (that is, the nuclide had captured more 
neutron after reaction as Eq. (8)), remain after cooling with- 
out decay because they are stable (Fig. 7). 

On the other hand, in the case of the creation of NRMs with 
odd atomic number Z’ (ex. 4;Rh), target FPs (,,Ru) have an 
even atomic number Z and many stable isotopes (,,Ru: 7). 
That is, even if an isotope (,,Ru) captures a neutron, it cannot 
decay easily because it probably changes another stable iso- 
tope, such as in Fig. 10. Fig. 11 shows the time variation of 
the weight of isotopes of target ,,Ru in the case of irradiation 
by LWR. It turns out that '}°*Ru, which has a possibility of 
creation of !™Rh, is the only stable isotope of ,;Rh by (n, y) 
reaction and B~ decay, as 


102Ru(n, y) 42Ru Aea 103Rh(stable), (9) 


had not decreased, but increased. It suggests that not !°Ru, 
but '°'Ru had mainly captured many neutrons and more nu- 
clei of !'?Rh had been created by (n, y) reaction of 'OlRu, 
then decreased them such as 


10TRu(n, y) 47 Ru(stable). (10) 


Thus, it turns out that the creation rate of ,,Rh was not high 
(Fig. 4) because target ,,Ru had been wasting many neutron- 
s by other reactions not involved in the creation of ,;Rh. In 
addition, Fig. 12 shows the time variation of the weight of 
isotopes of ,,;Rh created from target ,,Ru in the case of irra- 
diation by LWR. It turns out that all created isotopes of ,;Rh 
except 'Rh had decayed, although they had been created 
(Fig. 10). Therefore, many neutrons must be captured before 
they reach a B~ -unstable isotope which finally undergoes B7 
decay. In this case, the creation rate of a 7,,NRM element 
is not directly affected by a single neutron capture process, 
but is affected by the number of neutron capture reactions to 
reach to the ~ -unstable isotope of the target ,FP element. 
Since the number of stable isotopes varies from element to 
element, the dispersion from the average trend is expected to 
be larger in the latter case. 


IV. CONCLUSION 


A burnup calculation has been performed in order to ana- 
lyze the process which aims at creating highly-valuable stable 
elements (i.e. secondary nuclear rare metals; NRMs) from 
fission products (FPs) by using a neutron capture reaction, 


$10311-7 


Atsunori Terashima and Masaki Ozawa 


followed by B~ decay. As a result, high rates of creation of 
z+ıNRM elements were obtained for platinum group metal- 
s (44Ru: 8.03%/year (FBR), 4,Pd: 20.6%/year (LWR)) and 
middle or heavy rare earth elements (ex. ¢,Gd: 23.4%/year 
(LWR), DY: 17.1%/year (FBR)). It was found that creation 
rates of NRMs with even atomic number Z’ (= Z + 1) were 
generally high, but those of NRMs with odd Z’ were compar- 
atively low. Having introduced average 1-group neutron cap- 
ture cross sections of target ¿FP elements, correlation with 
creation rates of 7,,NRM elements was analyzed and the fol- 
lowing rules were established: the creation rate of a NRM el- 
ement with even (odd) atomic number Z’ was approximately 
proportional to the 0.78th (0.63th) power of the average 1- 
group neutron capture cross section of the corresponding tar- 
get FP element with odd (even) atomic number Z. It was also 
found that the dispersion from the average trend was much 
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larger in the case of creation of NRMs with odd Z’ than in 
the case of creation of ones with even Z’. These differences 
could be attributed to the even-odd effect of average 1-group 
neutron capture cross sections and the number of stable iso- 
topes. Based on the methodologies established and results 
obtained in this work, a design optimization of a reactor core 
and condition of neutron irradiation is going to be performed 
as a future work in order to increase the creation rate of N- 
RMs. 
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